Theory of the 7r state in 3 He Josephson junctions 
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The flow of superfluid 3 He-B through a 65 x 65 array of nanometer size apertures has been measured 
recently by Backhaus et al. They find in the current -phase relation a new branch, so-called n state. 
We study two limiting cases which show that the n state arises from coupling of the phase degree 
of freedom to the spin-orbit rotation. The n state exists in a single large aperture, but is difficult to 
observe because of hysteresis. A better correspondence with experiments is obtained by assuming a 
thin wall, where the Josephson coupling between the two sides arises from a dense array of pin-holes. 



PACS: 67.57.Np 

The flow of superfluid 3 He-B through a single 
nanometer-size aperture was studied by Avenel and Varo- 
quaux some time ago Q. At temperatures near the su- 
perfluid transition temperature T c , the current-phase re- 
lation is sinusoidal 



J{4>) = J c sin0 



(1) 



as expected for a Josephson junction. Also according to 
expectation, they find that the sine form (Q) gradually 
becomes tilted when the temperature is lowered. More 
recently, Backhaus et al. studied a 65 x 65 array of small 
apertures [^J^]. They discovered a new behavior where 
the current-phase relation acquires a positive slope at 
phase differences <fi ~ 7r - This 7r state develops when the 
temperature is lowered to approximately 0.6T C . 

A few theoretical explanations for the ir state have 
been proposed pL||. In this letter we present a theory, 
that is based on the many-component form of the order 
parameter in 3 He. It differs from the previous sugges- 
tions because it contains no unjustified assumptions and 
an order-of-magnitude agreement with experiments is ob- 
tained without any adjustable parameters. 

Unusual current-phase relations occur also in other sys- 
tems. A 7T junction, where J c in (Q) is negative, can be 
induced by adding magnetic impurities to a tunneling 
barrier between two s-wave superconductors p|. Sim- 
ilar 7r shifts appear in nonmagnetic junctions between 
d-wave superconductors. In addition, current-phase re- 
lations with additional zeros (J((f>) = for <f> ^ or 
7r) can appear for special orientations of the anisotropic 
crystals J7]. The ir state in 3 He differs from these in sev- 
eral respects, most fundamentally because it arises from 
interplay of two soft modes of the order parameter, the 
phase <j) and the spin-orbit rotation. 

We present calculations in two limiting cases. In the 
case of a tunneling barrier, the existence of the 7r state 
can be demonstrated by analytic calculations. The pa- 
rameters of the tunneling model are estimated using the 
quasiclassical theory. In the case of a single aperture, 
the 7r state is obtained by numerical simulations using 
the Ginzburg-Landau theory of 3 He. 

Tunneling junction. — The simplest case to demon- 



strate the 7r state is to consider a planar wall through 
which the 3 He atoms can tunnel. The energy arising 
from tunneling between the left (L) and right (R) sides 
can be written as Iral 



Fj = -Re£[aA^ + 6(^A 



R 
fix 



(2) 



Here A^j is the 3x3 matrix order parameter where 
the first index fi refers to the orientation of the Cooper 
pair in spin space, and the latter index j in orbital 
space. The z axis is taken perpendicular to the tun- 
neling wall. Equation (||) is a simple generalization of 
Fj = —a Re(A L * A R ) , which describes the Josephson cou- 
pling of two s-wave superconductors with order parame- 
ters A L and A R ||. The mass current through the wall 
is given by J = (2m^,/K)dF t ] / 'd(f>, where </> = <fr L — cf> R is 
the phase difference and UI3 the mass of a 3 He atom. 

In the B phase of 3 He, the order parameter has the 
form A^j = Aexp(i(f>)R^j. Here A is the amplitude, 
exp(i</>) a phase factor, and R^j a 3 x 3 rotation ma- 
trix: ^Z^R^jR^k = Sjk- The rotation matrices can be 
parametrized by an axis n and an angle 9. Substituting 
into (0) gives (a, (3 > 0) 



Fi 



•E 



ctR^zR^z 



)] cos</>. 



(3) 



In deriving this expression from (||) one must pay atten- 
tion to the fact that the order parameter of the p-wave 
superfluid is strongly suppressed near a wall. As a con- 
sequence the parameters a and (5 in (0) are not simply 
related to the coefficients a and b in (H), but otherwise 
the dependence of Fj on the soft variables (j) and 
remains the same as obtained by the simple substitution 
above. 

Let us consider the case that the rotation matrices 
on the left and right sides are the same. This gives 
rise to the "zero-state" with the critical current J c = 
(2m,3/h)(a + 2/3) > 0. This state has lowest energy 
when j 0| < 7r/2 because it corresponds to the maxi- 
mum of the expression in square brackets in (|^). The 
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FIG. 1. The current-phase relationships and energies for 
the tunneling model. The left and right panels correspond 
to parallel and antiparallel n's far away from the junction, 
respectively. The directions near the junction are depicted 
by arrows. The curves correspond to different values of the 
gradient-energy parameter 7: ideal it state (7 = 0, solid line), 
no 7T state (7 = 00, dotted line), and an intermediate case 
(7 = 0.245 aj, dashed line). The parameters a = 0.2207 aj 
and (3 — 0.0347 aj are chosen to imitate the experiment || 
at T = 0.55T C . 



situation changes when ci exceeds tt/2. There one has 
to look for a minimum of the expression in the square 
brackets. This corresponds to the it state, which is illus- 
trated by the solid line in Fig. [j]a. The critical current 
J c in ([!]) is negative: J c = —(2m 3 /h)a if a > and 
J c = —(2m3/h)(2f3 — a) otherwise. 

In order to make the tunneling model realistic, we 
have to consider three additional contributions to the en- 
ergy. Firstly, there is the magnetic dipole-dipole energy 
F d = 83 d A 2 (i + cos6>) 2 @. In the bulk it fixes the ro- 
tation angle 9 equal to 9q — arccos(— j) = 104°. This 
remains valid also near the junction because both the 
Josephson energy (||) and the dipole-dipole energy can 
reach their minima simultaneously: the products of two 
rotation matrices appearing in the former are not limited 
by the fact that both matrices have a fixed rotation angle 
9q. Secondly, there is a surface energy that arises from 
coupling of the dipole-dipole energy to the suppression 
of the order parameter near walls pT| . It has the form 



F s = 6 4 (n-s) 4 -6 2 (n-s) 2 , 



(4) 



where s is the surface normal. The lowest surface energy 
is achieved when the rotation axis n is perpendicular to 
the wall, n = ±s, because 62 > 264 > 0. Thirdly, there is 
a gradient energy associated with spatial bending of the 
rotation axis n. It arises because in practice all tunnel 



junctions are of finite size, and other walls in the con- 
tainer favor a different orientation of n than may be the 
minimum of the Josephson energy. We model the gradi- 
ent energy by the simple quadratic forms 



L \2 



1W - 



(5) 



where 77 is the polar angle of n, i.e., cos 77 = h z . rj L 
and r] R denote the polar angles on both sides just at 
the junction, and we assume that the values 77^ and 77^, 
further away are either or tt. In the experimental case 
B the surface energy (^) is important in fixing 77^. and 
77^ , but otherwise its contribution is so small that we can 
neglect it in the following. 

The current-phase relations for the tunneling model, 
(H) and (H), are plotted in Fig. [l] It can be seen that a 
large value of the gradient energy parameter 7 suppresses 
the 7T state. Furthermore, we find two cases where the 
rotation axes n far from the junction are either parallel or 
antiparallel. The latter has smaller critical current [ J c = 
{2m^/K){a — |/3)] but a relatively a more pronounced 7T 
state. This "bi-stability" was theoretically discussed in 
Ref. |l2| and has recently been observed experimentally 

i- 

Evaluation of tunneling parameters. — The experiment 
has a square array of apertures of diameter D = 100 
nm with spacing S = 3 [im in a wall of thickness W — 50 
nm. In order to make the tunneling model to imitate 
the experiment, we estimate a and (3 by letting all the 
three lengths to approach zero but keeping their ratios 
unchanged jl3| . The calculation for such "pin-holes" |l4|] 
is relatively simple once the self-consistent solution for 
the order parameter near a wall is known j^J. We as- 
sume diffuse scattering of quasiparticles at surfaces. The 
tunneling form (|^) is reproduced in these calculations at 
temperatures T > 0.5T C , and values of a and j3 can be 
extracted. The parameter 7 can be estimated using the 
bending energy of the B phase fl6f and assuming a simple 
form 77 (r) —7700 cx r , where r is the radius from the cen- 
ter of the aperture array, and this expression is cut off at 
the radius of the array. In agreement with experiments, 
we find that the 7T state appears at low temperatures be- 
cause a and /3 oc (1 — T/T c ) 2 have stronger temperature 
dependence than 7 oc 1 — T/T c , see Fig. g. Moreover, the 
parameters a, /3, and 7 agree within one order of mag- 
nitude to those that give an approximate best fit to the 
experiments (see caption of Fig. ||). This fit reproduces 
also the absolute magnitude of the critical current, and 
the same values of the parameters are used for cases of 
both parallel and antiparallel n's. 

The tunneling model can be improved trivially by ex- 
tending the pin-hole calculation to the whole temperature 
range < T < T c . A more ambitious project for the fu- 
ture would be the self-consistent calculation for aperture 
sizes on the order of the coherence length £o- In both 



cases the resulting Josephson energy Fj ( 
no more be of the simple form (|^). 



R R ) 



will 



2 




FIG. 2. The current-phase relationships for parallel (a) 
and antiparallel (b) n's far away from the junction. The differ- 
ent curves correspond to temperatures from 0.45T C to 0.85T C 
with intervals of 0.05T C . The parameters a and f3 are calcu- 
lated with the pin-hole model and 7 estimated as explained 
in the text. However, in order to get better correspondence 
with experiments, we have multiplied the estimated values by 
factors 7.5, 1.3, and 0.15, respectively. 
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FIG. 3. The order parameter in the tt state along the axis 
z of the aperture. The figure corresponds to the phase differ- 
ence (f> — Ti, i.e., the order parameters far left and far right 
differ by factor -1. The broken symmetry allows nonzero A yz 
and Ayz which have a long tail in the bulk. These compo- 




nents vanish in the zero branch. 
W = 6£gl and D = 10£gl. 



The wall is shown as shaded, 



FIG. 4. Theoretical phase diagram for the the tt state in a 
single aperture. The current and the energy as a function of 
phase difference <j> are shown in two cases. D is the diameter 
of the aperture and W the wall thickness. The tt branch is 
found in regions (b), (c), and (d), where it is locally stable 
at a fixed phase difference <j> ~ tt. In regions (c) and (d) 
it is also a local minimum of energy with respect to (f> at 
4> — tt. In regions (b) and (c) the tt state is the absolute 
minimum energy state at cj> = tt. The parameters of one 
aperture of the experimental aperture array [^j are shown by 
dashed line, and the observation of the tt state is marked by 
a cross. The temperature dependent GL coherence length 
is defined by £gl = Rwf/vI0A(T), where vf is the Fermi 
velocity and A(T) the BCS energy gap. 



Single aperture. — The limit opposite to the tunneling 
barrier is a single large aperture. There the major task 
is to calculate the order parameter self-consistently. We 
have done this using the Ginzburg-Landau (GL) theory of 
3 He. The differential equations were solved numerically 
on a grid in and around the aperture. Our calculations 
are more general than the previous ones |[2]|l7]] because 
we use a full a three-dimensional grid. Vanishing A^j 
was assumed at surfaces. 

The order parameter of the tt state is shown in Fig. 
|3|. It is plotted along the axis of a circularly symmetric 
aperture. For simplicity, we have normalized the order 
parameter to unit matrix in the bulk, A fl j{z = ±00) = 
exp(±i0/2)(5 M:) (assuming the case of parallel n's). This 
is possible because for aperture sizes on the order of the 
GL coherence length £gl the dipole-dipole energy can be 
neglected. The characteristic property of the tt state is 



the components A yz and A zy . These are the dominant 
components in the orifice, and they decay slowly towards 
the bulk. They imply broken symmetry: the symmetry 
group of the tt branch is m'm2' compared to 3™ of the 
zero branch. (Here prime denotes time-inversion.) This 
sets rather strong requirement for the calculation because 
the circular symmetry of the aperture cannot be used to 
simplify the computation. 

The current-phase relations are summarized in Fig. 
||. We see that the occurrence of the tt branch depends 
sensitively on the diameter of the aperture, whereas the 
wall thickness is less important. For small apertures no tt 
branch is found. When D exceeds approximately 5£gl, 
the tt branch appears. In the region (b) the current-phase 
relation has negative slope. Such a state can be stabilized 
if the left and right sides are connected, like in a torus 



3 



geometry. Increasing the diameter, the current-phase 
relation gets a positive slope in region (c). This state is 
stable also in a piston-driven flow channel. The tt state 
is also the absolute energy minimum at <f> = tt in region 
(c). In region (d) the tt state continues to exist but it 
has higher energy than the zero branch. The calculation 
assumes the idealized case of flow between two infinite 
bulk fluids. Any additional hydrodynamic inductance 
shifts upwards the border between regions (c) and (d) as 
it increases the energy of current carrying states. 

Although the tt branch constitutes the absolute energy 
minimum, it may be difficult to find it experimentally 
in a single aperture. The reason is that whenever it is 
locally stable, there always exists a locally stable zero 
state at the same 4>. Because the order parameters of 
the tt and zero states differ considerably, it may be that 
the phase slips only take place between two branches of 
the zero state without ever finding the way to the lower 
energy tt state. This is what we find in the numerical cal- 
culations, where the tt state was found only if the initial 
A^j was chosen close enough to the converged solution. 
We remind, though, that our calculations are not meant 
to simulate the correct dynamics of the phase slip. 

The dimensions of one aperture in the array at Berke- 
ley are marked on figure |J This is clearly in the re- 
gion where no tt state is found. Although the Ginzburg- 
Landau calculation is accurate only at temperatures near 
T c , it is unlikely that the tt state could be stabilized in 
a more accurate calculation at lower temperatures. This 
statement is based on experience gained in previous low- 
temperature calculations pl|. Thus we conclude that 
the appearance of the tt state in the Berkeley experiment 
essentially depends on the presence of many apertures. 

We have also done two-dimensional calculations that 
simulate the flow through a long narrow slit. The tt state 
is found and its properties are qualitatively similar to 
those in a circular aperture fll9| . In particular, the tran- 
sitions from the zero branch to the tt branch seem to be 
absent. This is consistent with the fact that no tt branch 
was found in the experiments by Avenel and Varoquaux 

The tt state can be interpreted so that a half- quantum 
vortex has crossed the orifice. There are no free half 
quantum vortices in superfluid 3 He-B, but the double- 
core vortex ]2Q] can be interpreted as a bound pair of two 
half-quantum vortices ]21j| . Indeed, the order parameter 
in Fig. is very similar to that in the double-core vortex 
on the axis going between the two cores pc| . 

Conclusion. — The tt state was found to occur in both 
the limits investigated above. Its mechanism is the same 
in both cases: a lower coupling energy is achieved by 
producing a spin-orbit rotation that heals slowly in the 
bulk liquid. The Berkeley experiment O is somewhat an 
intermediate case to the two limits studied, which makes 
it evident that the tt state also there arises from the same 
mechanism. 



The theory above provides several predictions that can 
be tested experimentally. For example, the tt state de- 
pends on the linear dimension L of the aperture array 
because a/7, (3/j oc L. An external magnetic field fixes 
the surface orientation of n, and thus can be used to 
suppress the tt state. 
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